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Using Gutzwiller-projected wave functions, I estimate the ground-state energy of the t-J model 
for several variational states relevant for high-temperature cuprate superconductors. The results 
indicate antiferromagnetism and phase separation at low doping both in the superconducting state 
and in the staggered-flux normal state proposed for the vortex cores. While phase separation in the 
underdoped superconducting state may be relevant for the stripe formation mechanism, the results 
for the normal state suggest that similar charge inhomogeneities may also appear in vortex cores up 
to relatively high doping values. 



After many years of research, the issues of antifer- 
romagnetism and of phase separation in weakly-doped 
high-T c cuprates are far from being settled. While it 
is generally accepted that short-range antiferromagnetic 
(AF) correlations are crucial for the superconductivity 
in cuprates, it is not clear whether the long-range anti- 
ferromagnetism at very low doping should be considered 
as a competing order on equal footing with the super- 
conductivity or simply as a minor side effect. The latter 
point of view is implicitly assumed in the resonating- 
valence-bond (RVB) scenario of high-temperature super- 
conductivity [1]. In terms of ground-state properties, the 
RVB approach may be recast into the language of vari- 
ational Gutzwiller-projected (GP) wave functions [2-4]. 
Recently a considerable progress has been reported in de- 
scribing properties of high-temperature superconductors 
with the help of GP wave functions for the t-J (or Hub- 
bard) model [5]. Most studies of the GP wave functions 
neglect the long-range AF ordering at low doping, and 
the resulting phase diagram contains the superconduct- 
ing phase starting from zero doping. However, it is very 
easy to take AF order into account by explicitly including 
it in the variational wave function. Within this approach, 
the GP wave function for the t-J model is known to be 
energetically unstable with respect to the AF order below 
the level of doping about 10% (at J/t = 0.3) [6]. 

The phase-separation issue is a much more delicate 
subject than antiferromagnetism: it is not decided even 
at the level of the t-J model. While the phase separa- 
tion at large J/t is well established, different studies do 
not agree about whether the phase separation occurs in 
the physically relevant parameter range (at J/t ~ 0.3) 
[7-10]. 

In this work I refine the numerical results of Ref . [6] on 
the ground-state energy of the superconducting state in 
presence of AF order and perform a similar analysis for 
the staggered-flux state recently proposed to describe the 
normal state in the vortex cores [11,12]. For both normal 
and superconducting states we find antiferromagnetism 
and phase separation at low doping. The phase separa- 
tion follows from the upward convexity of the ground- 
state energy as a function of doping. Within our approx- 



imation, the phase separation persists to higher dopings 
than antiferromagnetism. Consequently, the coexistence 
of AF order and superconductivity is not realized as a 
homogeneous state (at least, for the considered dimen- 
sionless parameter t/J — 3). Instead, at dopings lower 
than the phase-separation point x scp , it is energetically 
favorable to split into two phases: the undoped antifer- 
romagnet (with the long-range AF order but without su- 
perconductivity) and the superconductor with the doping 
x scp (without the AF order). The relative areas of the 
two domains are fixed by the average doping, and the ac- 
tual shape of the domains is determined by the Coulomb 
interaction at large distances (neglected at the level of 
the t-J model). This is the phase-separation scenario 
of the stripe formation [13-15]. Some numerical studies 
suggest an alternative point of view that stripes appear 
already in the t-J model without any long-range interac- 
tions [16]. The results reported in this paper do not have 
any implications on the latter scenario, since I do not con- 
sider here any incommensurate spin-charge-density-wave 
states. A priori it is possible that stripe states with ener- 
gies lower than the states studied here exist and may also 
be constructed variationally by Gutzwiller projection. 

Regarding the effect of the long-range Coulomb inter- 
action, I only remark that already nearest-neighbor re- 
pulsion may act to reduce the phase separation. With in- 
creasing nearest-neighbor repulsion, the phase-separation 
region shrinks to smaller doping values and practically 
disappears at the repulsion strength about 3-4 J — ap- 
proximately equal to the repulsion strength required to 
suppress superconductivity. [Note that nearest-neighbor 
repulsion may also favor stripe-like states not considered 
here.] 

Finally, I comment on a similar phase-separation fea- 
ture in the staggered-flux normal state proposed recently 
to describe the normal state in the vortex core [11,12]. 
The phase separation in the normal state appears more 
prominent and extends to higher doping values than in 
the superconducting state. This suggests that charge in- 
homogeneities (an analogue of stripes) may also appear 
in the vortex core. However, as estimated in Ref. [12], 
the core size should be of order of several lattice spac- 
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ings, and large gradients of the order parameter play an 
important role in determining the structure of the core, 
together with the long-range Coulomb interaction. Thus 
the problem of the vortex core structure, even at the 
level of optimizing the Gutzwiller-projected wave func- 
tion, becomes a very complicated one, and our present 
analysis of uniform states is not sufficient for solving it. 
We note here that insulating and antiferromagnetic re- 
gions in the vortex cores have also been predicted in other 
approaches: in the SO(5) theory [17] and in the unre- 
stricted mean-field analysis [18]. 

In the rest of the paper, I present the details of our 
numeric analysis of the Gutzwiller-projected wave func- 
tions. 

The variational wave functions studied in this work are 
chosen to minimize the energy of the t-J Hamiltonian 



H = P G 



Pg , (1) 



where the sum is taken over the pairs of nearest- 
neighboring sites i and j on the two-dimensional square 
lattice, rii denotes the hole density at a given site. Fol- 
lowing the usual Gutzwiller-projection approach, we con- 
sider the variational wave functions of the form V&gp = 
Pg^o- Both here and in the Hamiltonian (1) Pg denotes 
the "double" projection: it projects out components with 
doubly occupied sites (the usual Gutzwiller projection) 
and further it also fixes the total number of particles to 
the required value, Wo is the ground state of the auxiliary 
("mean-field") BCS Hamiltonian 

Hbcs = J2 (-XtfCfaCja + A ij(4Al ~ c !| c Jt) + h - c -) 

ij 

+ J2(-iyhaf, (2) 



are estimated to be of the order of the Monte Carlo sta- 
tistical errors (about J) [in the staggered-flux state 
the finite-size corrections are somewhat bigger because of 
the discretization of the Fermi pockets] . The parameter 
of the t-J model is taken to be t/J = 3 throughout the 
paper. 

To check for an instability with respect to the AF or- 
dering, we further minimize the variational energy as a 
function of h while keeping A/x and [i/x constant (in 
principle, the minimization should be performed by vary- 
ing all the three parameters simultaneously; we can how- 
ever check that, in the vicinity of the energy minimum, 
the errors from our simplified minimization procedure are 
negligible) . Instead of characterizing the AF state by the 
fictitious field h, we use a physically significant quantity, 
the staggered magnetization. In Fig. 1 we plot the stag- 
gered magnetization in the optimal wave function as a 
function of doping x. We find the instability towards 
the AF ordering below xaf =0.11. So far the compu- 
tations repeat those in Ref. [6]. Thus obtained values 
of staggered magnetization should only be considered as 
an indication of an instability, and not as good numer- 
ical estimates of the actual staggered magnetization in 
the ground state: for example, in the undoped case, the 
variational estimate is much higher than the known exact 
value [19]. 
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where of = cLc^ — cjiCjj, is the z-magnetization at the 
site i. As usual, we take the nearest- neighbor d-wave 
ansatz: Xij = X on nearest-neighbor links, xu — A* is the 



on-site chemical potential, Aj 



±A, with ± for ver- 



tical and horizontal nearest-neighbor links respectively. 
Within this variational ansatz, the wave function depends 
on the three dimensionless parameters: A/x, /u/x> and 
h/x- The parameter h represents an artificial staggered 
magnetic field acting on spins to produce the long-range 
AF order. At zero value of h, the minimization has been 
performed by many authors [3-5,12]. We further adopt 
the numbers reported in our earlier publication Ref. [12]. 
As in that work, we compute the ground-state energy by 
the variational Monte Carlo method (see e.g. Ref. [3] for 
details of the method), on the square lattice 22x22 with 
the boundary conditions periodic along one direction and 
antiperiodic along the other direction. For this system 
size, the finite-size corrections to the ground-state energy 



FIG. 1. Staggered magnetization M = ((— l)Vf ) as a func- 
tion of doping at t/J = 3 for the superconducting (circles) and 
staggered-flux (squares) states. Vertical error bars are of the 
order of the symbol size. Error bars of the phase transition 
points are about 0.01 (shown). Our results for the supercon- 
ducting state agree with those in Fig. 2 of Ref. [6] . The arrow 
shows the exact value (M — 0.60) of the staggered magneti- 
zation in the undoped Heisenberg antiferromagnet. 

If we analyze the resulting ground-state energy as a 
function of doping, we observe that it is upward convex 
at low doping. This implies that at those dopings phase 
separation occurs. In Fig. 2 we plot the ground-state en- 
ergy with the subtracted linear part E — /i S cp^ as a func- 
tion of doping. The slope /j, scp of the limiting tangent 
line gives the chemical potential at the phase separation 
point. The point of contact determines the critical dop- 
ing x sep below which the phase separation occurs. At 



t/J = 3, wc estimate x. 



0.13 and 



/^scp 



-5.65J 
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(this value of x scp agrees with the estimates of Ref. [8] ) . 

Note that for the t-J model with t/J = 3 we obtain 
x scp > xaf, i-e. antiferromagnetism does not mix with 
superconductivity in a uniform phase; instead, phase 
separation occurs between the undopcd antiferromagnet 
and the superconductor at doping x scp without antifer- 
romagnetism. However, this inequality appears to be 
non-universal, and could, in principle, be reversed by 
adding other terms like next-nearest-neighbor hopping or 
nearest-neighbor repulsion (see the discussion of nearest- 
neighbor repulsion below). 

We may further perform the same steps as above for 
the "normal" staggered-flux (SF) state proposed to de- 
scribe vortex cores in the mixed state [11,12]. We first 
test the SF state for the AF instability and find that, 
similarly to the superconducting state, the SF state fa- 
vors AF order at very low doping. We find that the 
corresponding staggered magnetization (also plotted in 
Fig. 1) is smaller than in the superconducting state at 
the same doping. The critical value of doping x^^ is 
also smaller than in the superconducting state (we esti- 
mate x^p" 1 — 0.08). Thus antiferromagnetism in vortex 
cores appears at the first glance more fragile than in the 
bulk of the superconductor. 




FIG. 2. Energies of different GP wave functions with linear 
part subtracted (E — fi scp x) as a function of doping x. The 
five different variational wave functions are compared (opti- 
mized in energy within each class): superconducting without 
antiferromagnetism (solid circles, solid line), superconducting 
with antiferromagnetism (empty circles, dashed line), stag- 
gered-flux without antiferromagnetism (solid squares, solid 
line), staggered-flux with antiferromagnetism (empty squares, 
dashed line), and zero-flux (projected Fermi sea, empty dia- 
monds, dotted line). Left panel (a) shows energies in the 
absence of NN repulsion (V = 0), with /i scp = —5.65. Right 
panel (b) corresponds to V = 4J, with /i scp = —4.85 [NN 
repulsion is included at the perturbation level via (3)]. The 
energies are in the units of J per lattice site. Error bars are 
smaller than the symbol size. For comparison, the arrow in 
panel (a) shows the exact ground-state energy in the undoped 
case (E = -0.669 J) [19]. 

However the phase-separation effect in the SF state is 
more pronounced than in the superconducting state: we 
estimate ^ = 0.21 — far in the "overdoped" region 



of the phase diagram. We summarize our results on the 
ground-state energies of different states in Fig. 2 (the 
difference between the energies of the superconducting 
and SF states was reported previously in Ref. [12] as the 
condensation energy for the superconducting state) . Our 
results on the phase separation in the SF state indicate 
that charge inhomogencitics (similar to stripes in under- 
doped cuprates) are likely to form in the vortex cores 
up to rather high doping values. As in the case with 
stripes, the actual structure of inhomogeneities should 
be determined by the long-range Coulomb interactions. 
Such inhomogeneities may play a role in producing the 
"subgap state" features observed in the density of states 
in the STM experiments [20] . 

As a simplified version of long-range Coulomb inter- 
action, we consider the effect of nearest-neighbor (NN) 
repulsion on the phase separation. The NN repulsion is 
included as the additional term in the Hamiltonian (1), 
H\r = V riiUj , where the sum is taken over all NN 
pairs of sites [n, denotes the hole density, as before; V is 
the repulsion strength in addition to the attraction — J/4 
already present in the Hamiltonian (1)]. If wc neglect in- 
teraction between holes, Hy might be thought to scale 
as x 2 with doping x. Then a small repulsion energy V 
might be already sufficient to change the convexity of the 
total energy as a function of doping and hence to sup- 
press phase separation. In reality, in the superconducting 
antifcrromagnetic state at small doping, the attraction 
between holes is quite strong, and scaling of Hy with 
doping is closer to x than to x 2 . In Fig. 3a, b we plot the 
nearest- neighbor correlation function (niUj) normalized 
by x and by x 2 , respectively. 
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FIG. 3. Nearest-neighbor correlation of the hole density 
(niiij) as a function of doping for the same wave functions 
as in Fig. 2. Panel (a) shows nearest-neighbor correlation 
divided by the average hole density x, and in panel (b) the 
same correlation is divided by x 2 . 

To estimate the effect of NN repulsion, we approximate 
the ground-state energy with repulsion as 

E(x)=E (x) + 2V(n i n j ), (3) 

where E (x) is the variational ground-state energy with- 
out repulsion, and the average is taken over the varia- 
tional ground state without repulsion. In this approxi- 
mation, we can estimate the repulsion strength necessary 
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to reverse the convexity of the ground-state energy E{x). 
With increasing the NN repulsion strength V, the upward 
convexity of the ground-state energy diminishes, and the 
phase-separation point a; sep shifts towards smaller dop- 
ing (both in the superconducting and the SF states). At 
the same time, the superconductivity gets also strongly 
reduced (at least in its present nearest-neighbor d-wave 
ansatz), with the superconducting transition point shift- 
ing to smaller doping. The antifcrromagnetism at finite 
doping is also reduced, but not so strongly as the super- 
conductivity. 

With increasing the NN repulsion strength V, the 
phase separation eventually becomes undetectable within 
our numerical errors at V ~ 3.. 4 J. At the same time, at 
the repulsion strength of about 4J, our superconduct- 
ing wave function looses the energy competition to the 
staggered- flux one: the price of nearest-neighbor holes 
overweighs the energy gain from the superconductiv- 
ity. It is possible however that at such strong NN re- 
pulsion another superconducting ansatz (e.g., involving 
next-nearest-neighbor pairing) would be more energeti- 
cally favorable (search for such states goes beyond the 
scope of the present paper). We illustrate the effect of 
NN repulsion in Fig. 2b where the marginal situation is 
shown: the superconducting state is nearly equal in en- 
ergy to the SF state at V — 4 J [with NN repulsion taken 
into account only to the lowest order in the perturba- 
tion theory (3)]. Note that our numerical results show 
a rather flat doping dependence of the ground-state en- 
ergy at very low doping, and this leaves unresolved the 
issue of whether phase separation at infinitesimally small 
doping may survive arbitrarily strong NN repulsion. Of 
course, the above treatment of the hole interaction is only 
a rough estimate. As suggested in Ref. [21], the varia- 
tional wave function may be further improved by includ- 
ing additional Jastrow-type factors. Such a modifica- 
tion of the wave function apparently affects the hole-hole 
correlations [21], and may therefore be important for a 
proper assessment of the effects of the hole-hole interac- 
tion. 

The variational study reported in this paper should be 
taken with care when applied to actual cuprate super- 
conductors: the phase-separation effect in the supercon- 
ducting state is weak, and many additional factors may 
change the picture outlined here. However, I believe that 
some qualitative conclusions may withstand small pertur- 
bations of the model: (1) staggered-flux phase is more 
disposed to phase separation than the superconducting 
state. This may produce charge inhomogeneities inside 
the vortex cores; (2) antifcrromagnctic ordering at low 
doping enhances phase separation (this may be a general 
property of phase separation in doped Mott insulators, 
see, e.g., Ref. [22] and references therein); (3) both in the 
staggered-flux state and in the superconducting state, the 



long-range antifcrromagnctic ordering occurs only at low 
doping (below ~0.1). Thus static AF order in vortex 
cores [23], if confirmed, possibly indicates regions with 
reduced hole concentration. 
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